Compressibility of Interacting Electrons in Bilayer Graphene 
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Using the renormalized-ring-diagram approximation, we study the compressibility ol the interact- 
ing electrons in bilayer graphene. The chemical potential and the compressibility of the electrons 
can be significantly altered by an energy gap (tunable by external gate voltages) between the va- 
lence and conduction bands. For zero gap and a typical finite gap in the experiments, we show both 
systems are stable. 
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Bilayer graphene has attracted considerable atten- 
tion because of its promised application in electronic 
devices^ - — In contrast to the Dirac fermions in mono- 
layer graphene, the energy bands of the free electrons 
in bilayer graphene are hyperbolic and gapless between 
the valence and conduction bands. 5 Most importantly, 
an energy gap opening between the valence and conduc- 
tion bands can be generated and controlled by external 
gate voltages. At low carrier doping, because the energy- 
momentum dispersion relation around the Fermi level 
is relatively flat, the Coulomb effect is expected to be 
significant in the interacting electron system of bilayer 
graphene. The Coulomb effect has been studied by a 
number of theoretical works using the Hartree-Fock (HF) 
and random-phase approximations.— ~— One of the ther- 
modynamic quantities directly reflecting the Coulomb ef- 
fect is the electronic compressibility. Recently, this quan- 
tity has been measured by experiments.— ~— For timely 
coordinating with the experimental measurements, it is 
necessary to theoretically study the combined effect of 
Coulomb interactions and the gap opening in the com- 
pressibility with a more realistic model. 



In this work, using the renormalized-ring-diagram ap- 
proximation (RRDA) ) 20 i 21 we study the Coulomb effect 
in the compressibility of the interacting electrons in bi- 
layer graphene. The chemical potential and its deriva- 
tive with respect to the carrier density are calculated 
for systems of zero gap and a typical gap in the exist- 
ing experiments. Though a finite-gapped system is more 
perturbable by the Coulomb interaction than the zero- 
gapped one at low carrier doping and low temperature, 
we show both systems studied here are stable. 



The atomic structure of bilayer graphene is shown in 
Fig. 1. The two sublattices in each layer are denoted by 
A (white) and B (black) atoms, respectively. The inter- 
layer distance is c = 3.34 A w 1.4a where a ~ 2.4 A is 
the lattice constant of monolayer graphene (the distance 
between two nearest A atoms). The energy of electron 
hopping between the nearest- neighbor (nn) carbon atoms 
in each layer is t f=s 2.82 eV,— while the interlayer nn hop- 
ping is h « 0.39 eVM 



The Hamiltonian describing the electrons is given by 
h = Yj e i n i ~ X! kA^io + SriiVijSnj (1) 

j ijcr ij 

where ej = ±A for top (back) layer, rij is the electron 

density at site j, c| ff creates an electron at site % with spin 
cr, 5rij = rij—n with n the average occupation number of 
electrons per site (which is also the charge number of the 
neutralizing background), My is the Coulomb interaction 
between electrons at sites i and j. Here A is the energy 
gap parameter, which is a consequence of the potential 
difference between the top and back gates (attached to 
the top and back layers, respectively). The model is re- 
stricted to nn hopping within the same layer and between 
the adjacent sites on top and back layers as shown in Fig. 
1. The Coulomb interaction is given by Vij — e 2 /ery for 
^ with e ~ 4 the dielectric constant of high fre- 
quency limit of the system; r,*j here is the distance be- 
tween two electrons at sites i and j in the bilayer system. 
For onsite Coulomb interaction U — Vjj , we take U in our 
calculation as U — 2vab where AB denotes the nn sites 
in the same layer. The dimensionless Coulomb coupling 
constant is given by e 2 /eat ps 0.5. As described by Eq. 
(fl"j), we here consider only the charge-charge interactions. 
The onsite interaction may lead to the antiferromagnetic 
correlation. We here neglect it since it is not important 
in graphene. 

By defining the doped electron concentration S per car- 
bon atom, we have n = 1 + 6. Under the transformation 
5 — > —5 and Cj >(7 <H> ±cj a for electrons at A (B) site, H is 
unchanged. Therefore our system satisfies the particle- 
hole symmetry. Furthermore, K = H — /j,(N — Nq) (with 
N the operator of the total number of electrons and Nq 
the total number of lattice sites and thereby N — No as 
the total number of doped electrons) is also unchanged 
under the above transformation and fi — > —fi. Thus, the 
chemical potential [i is an odd function of 5. 

To proceed, we begin with the Green's function G de- 
fined as 



(2) 



G(i,j,r-r') = -{T T C icr (T)Cl(r')) 

where C) a = (.c\ l2a ,c\ lja ,c\ 2ja ,c\ 2ja ) with c\ 
ating an electron of spin a at site A (B) of /zth (= 1,2 



ere- 




FIG. 1: Structure of Bernal stacking bilayer graphene. The 
energies of intra- and interlayer electron hopping between the 
A (white) and B (black) atoms are given by t and ti, respec- 
tively. 



respectively for top and back) layer of jth unit cell. In 
the momentum-frequency space, G is expressed as 



G(k, iu>i) = - 6fe - S(fc, 



with 




(3) 



(4) 



where uji = (21 + 1)tvT is the fermionic Matsubara fre- 
quency with I as integer number and T the temperature, 
£fc = 4o[l+exp(— ik x ) +exp(— ik y )], /x is the chemical po- 
tential, and E(fc, iu>e) is the self-energy. Under RRDA, £ 
is geven by the functional derivative of the 'free energy' 
functional <I> with respective to G shown diagrammati- 
cally in Fig. 2, 



S = 5&/SG. 



(5) 



In Fig. 2, each bubble is composed by two renormalized 
Green's function G's. The last terms in the expression 
of $ and S stem from the additional term in writing 
the product of two density operators appeared in the in- 
teraction term of Eq. ([T]) in normal order of electron 



SijCj^Cja. Note 



Operators: C^C^C^Cjc = C^C^CjaCia 

that Eq. J5) is a 4 x 4 matrix equation. In terms of G, 
the elements of £ are expressed as 



^fiu(k, iu>t) = -— G fll/ (k - q, iut - ^ m )W M „(g, iv„ 



where M = Nq/4 is total number of unit cells in the 
system, v m is the bosonic Matsubara frequency, and 



NU/2 




oo 



s = 



U/2 



FIG. 2: (color online) Diagrammatic expressions for 'free en- 
ergy' functional $ and self-energy E. The line represents the 
Green's function and the wavy line is the Coulomb interac- 
tion. The last terms stem from the additional term in writing 
the product of two density operators in normal order of elec- 
tronic operators. U is the onsite Coulomb interaction. 



W^, v (q,ip m ) is an effective interaction. The matrix form 
of W is given by 

W(q,iv m ) = [1 - v(q)x(q,iv m )r 1 v(q) 



with 



2T 
~M 



J~] Giiu [k, ioj£)G Ufl (k - q,ioj£ - iv m ) 



k-.e 



and v(q) is the Fourier component (matrix form) of the 
Coulomb interaction. In the present lattice model, v(q) 
cannot be expressed explicitely. A similar calculation of 
v(q) for the single layer graphene is given in Ref. l2ll . The 
chemical potential /i is determined by 



2T 

= — ^2 TrG(fc, iu) t ) exp(zw f ?7), 



Ni 



(6) 



k-C 



where rj is an infinitesimal small positive constant. The 
Green's function so determined satisfies the microscopic 
conservation laws<21 

The behavior of the chemical potential is closely re- 
lated to the compressibility k of the doped carrier system. 
It is defined as 



1 dn 



(7) 



Its inverse (d/i/dn)T can be calculated by numerical 
derivative of /x determined by Eq. ^ with respect to n. 
Clearly, when (dfi/dn)T goes to zero, K becomes infinity, 
implying an inhomogeneity tendency in the system. 

Figure 3 shows the chemical potential /x as a function 
of 6 for the systems of A = (left panel) and A/t = 0.02 
(right panel) at various temperatures. The results by 
RRDA are compared with that by the self-consistent 
Hartree-Fock approximation (SCHF) (which is obtained 
from $ by neglecting all other ring terms except the first 
one) for the zero-gapped system. The chemical poten- 
tial [i for zero gap case goes to zero as 6 — > because 




<?(10" 3 ) <?(10- 3 ) S(\Q- 3 ) 



FIG. 3: (color online) Chemical potential fi as function of 
electron doping concentration 8 at various temperatures cal- 
culated by SCHF and RRDA for zero-gapped systems (left 
panel) and by RRDA for systems of A/t — 0.02 (right panel). 



FIG. 4: (color online) (dfi/dn)T as function of doping con- 
centration S at various temperatures obtained by SCHF for 
A = 0. The inset shows the result at T/t = 10~ 2 for a wide 
range of doping. 



of the particle-hole symmetry. The overall behaviors of 
/i given by RRDA for system of A = at various tem- 
peratures are smooth and seem not much different from 
each other. Physically, at low doping, the electrons in 
the valence band can be thermally excited to high levels 
in the conduction band without significant change in the 
chemical potential because the density of states is sym- 
metric about the touch point of the valence and conduc- 
tion bands. The same feature is seen by SCHF. However, 
in contrast to the smooth behavior given by RRDA, the 
curve corresponding to lower temperature given by SCHF 
is twisty in a region close to zero doping. This behavior 
will result in a singularity in the compressibility. On the 
other hand, for the gapped system, the chemical poten- 
tials are indistinguishable only at very low temperatures. 
At low temperature, /i is finite at S = 0. For free electron 
system, fj, approaches A as S — > + and T — > 0. The limit 
here is less than A because of inter-electronic Coulomb 
interactions. Since the limits S — > ± are different, there 
is a discontinuity in fj, at S = from hole side to electron 
side. With increasing the temperature, fx is suppressed 
due to thermal excitations between the valence and con- 
duction bands. The thermal effect is significant at lower 
carrier doping since where the Fermi energy is lower. At 
high temperature, fi at 5 = + is bound from below, 
fi > 0, because of the particle-hole symmetry. 

For illustrating the RRDA result, we firstly show in 
Fig. 4 the one by SCHF for (d/j,/dn)T as a function of 
S for A = at various temperatures. At low temper- 
ature, there is a peak at 5 = 0. This is because the 
chemical potential at finite doping is suppressed by the 
Coulomb interactions but not changed at S = due to 
the particle-hole symmetry. At low but finite doping 
close to S = 0, there is a sharp decrease in (dfi/dn)T 



as T — > 0. This stems from the twisty behavior in /i. At 
T = 10~ 4 , (dfx/dn)T becomes negative for electron/hole 
concentration within a small region. At T = 0, the re- 
sult (dfi/dn)T < at low doping has been predicted 
by Kusminskiy et al. using pcrturbative HF approach 
on a continuum model— The compressibility is singu- 
lar at the doping where (d[i/dn)T — 0. This is a strong 
Coulomb effect at low doping where the Fermi level of the 
electrons is low. This is similar to the prediction of the 
Wigner crystallization in the two-dimension electron gas 
at low density. The inset shows the result at T/t = 1CP 2 
for a wide range of doping. All curves approach the same 
behavior at large doping as the one at T/t = 10~ 2 . 

The differences between Ref. [25| and the SCHF cal- 
culation here are in following aspects. (1) Ref. [H uses 
the perturbative HF in which the Green's function not 
renormalized as done in SCHF. (2) The lattice model 
here is different from the continuum model adopted in 
Ref. [25l In the continuum model, the intralayer and 
interlayer interactions are given by v(q) — 2i:e 2 /q and 
2ire 2 exp(— cq)/q, respectively, with a cutoff q c . Because 
of this cutoff, the interactions in real space are not the 
Coulomb form e 2 /rij that is used here. The compressibil- 
ity k should sensitively depend on the interactions. (3) 
By the lattice model, the particle-hole symmetry is satis- 
fied. On the other hand, there is no such a symmetry in 
the continuum model. (4) The present SCHF calculation 
is for finite temperature, different from Ref. that the 
calculation is performed for T = 0. 

In general, both HF and SCHF for the systems of 
long-range Coulomb interactions over count the strong 
coupling effects. More terms arising from the interac- 
tions should be included in the 'free energy' functional. 
By RRDA, all the ring terms that are the most long- 
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FIG. 5: (color online) (djj,/dn)T as function of doping con- 
centration S at various temperatures obtained by RRDA for 
A = 0. The inset shows the result at T/t = 10~ 2 for a wide 
range of doping. 



wavelength divergent terms are summed over in $. As 
a result, the effective interactions between electrons are 
weakened from the bare Coulomb one. Shown in Fig. 5 
is (dfi/dn)T obtained by RRDA as a function of S for 
the zero gap system at various temperatures. The re- 
sults given in Fig. 5 are quite different from that in Fig. 
4. At temperatures down to T/t = 5 x 10~ 4 , (d/j,/dn)T 
monotonically decreases with the electron/hole doping 
concentration. At T/t < 2 x 10~ 4 , the Coulomb suppres- 
sion shows up nearby the central peak similarly as in Fig. 
4. But the deviation from the monotonically decreasing 
behavior is much smaller than that given by SCHF. From 
the low temperature behaviors of (dfi/dn)?, we deduce 
that it is positive at T — and the system is stable. 

The results for a gapped system of typical A/t — 0.02 
as in experiments arc depicted in Fig. 6. By compar- 
ing to the zero gap case, the behavior of (d/j,/dn)T is 
delicate. At temperature higher than the gap, because 
there are significant excitations between the valence and 
conduction bands, both zero and finite-gapped systems 
are not different so much. However, at low temperature, 
(dfi/dn)? is suppressed much in a wide doping range 
around 6 = 0. This is because of the special conduction 
band structure of the gapped system. For A = 0, the 
band bottom is parabolic. At finite A, however, the lo- 
cation of the minimum of the conduction band (for elec- 
tron doping) is a ring surrounding the Dirac point. In 
a region encircling the ring, the energy band is nearly 
fiat. Therefore, at low doping, the energy of electrons is 
approximately dispersionless. For such a dispersionless 
band, the free electron system is more perturbable un- 
der the Coulomb interactions. Even though the quantity 
(d(x/dri)T is considerably suppressed, by carefully look- 
ing at its behavior at low T, it is positive in the limit 



FIG. 6: (color online) (d/j,/dn)T as function of doping con- 
centration S at various temperatures obtained by RRDA for 
A = 0.02 gapped system. 



T -> 0. 

The (dfi/dn)T < instability is not observed in recent 
experiments on the compressibility.— ~— In the experiment 
of Ref. |2j, the quantity (dfi/dn)T is observed as mono- 
tonically deceasing in both sides of the electron and hole 
doping with a central peak around the neutral point. In 
Ref. 3, the observation at zero magnetic field shows that 
the central peak in (d^i/dri)T is very sharp and the quan- 
tity then varies nearly flat within a wide range of elec- 
tron/hole doping. While in the experiment of Ref. El, 
the structure of central peak-dip-hump in (d^/dn)T as 
function of doping is observed (at zero magnetic field). 
Though the experiments do not quantitatively agree with 
each other, the common fact is (dfi/dn)? > 0. As for 
the dip or the flat region nearby the central peak in 
(d/i/dn)T, from the opinion of the present calculation, 
they may stem from the Coulomb suppression effect as 
studied here. However, the calculation given here is not 
in quantitatively agreement with the experimental re- 
sults. So far there is no theory can quantitatively explain 
the experiments. To solve the problem, besides the long 
range Coulomb interaction, other kinds of effects may 
need to be taken into account. 

In summary, we have investigated the compressibil- 
ity of the interacting electrons in bilayer graphene us- 
ing RRDA. For the zero gap system, with comparing to 
the result by SCHF, the compressibility given by RRDA 
is positive. The homogeneous system with zero gap is 
therefore stable. The system with a finite-gap is more 
perturbable by the Coulomb interaction than the one of 
zero gap. For the typical gapped system as investigated 
in experiments, however, we find that the system is stable 
as well. 
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